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This work describes the statistics for the occupation numbers of quantum levels in a large isolated 
quantum system, where all possible superpositions of eigenstates are allowed, provided all these 
superpositions have the same fixed energy. Such a condition is not equivalent to the conventional 
microcanonical condition, because the latter limits the participating eigenstates to a very narrow 
energy window. The statistics is obtained analytically for both the entire system and its small 
subsystem. In a significant departure from the Boltzmann-Gibbs statistics, the average occupation 
numbers of quantum states exhibit in the present case weak algebraic dependence on energy. In the 
macroscopic limit, this dependence is routinely accompanied by the condensation into the lowest 
energy quantum state. This work contains initial numerical tests of the above statistics for finite 
systems, and also reports the following numerical finding: When the basis states of large but finite 
random matrix Hamiltonians are expanded in terms of eigenstates, the participation of eigenstates 
in such an expansion obeys the newly obtained statistics. The above statistics might be observable 
in small quantum systems, but for the macroscopic systems, it rather reenforces doubts about self- 
sufficiency of non-relativistic quantum mechanics for justifying the Boltzmann-Gibbs equilibrium. 

PACS numbers; 



I. INTRODUCTION 

It is known empirically since the introduction of Quan- 
tum Hypothesis by Planck that thermal Boltzmann- 
Gibbs distribution works extremely well for quantum sys- 
tems. However, purely quantum derivation of this distri- 
bution is still not on a satisfactory ground. In a self- 
contained derivation one should be able to start from 
a large isolated system and then obtain the statistical 
distribution for a small subsystem. The conventional 
derivation of the Boltzmann-Gibbs distribution proceeds 
by postulating the micro-canonical condition. This con- 
dition has different status in classical and quantum me- 
chanics. In classical mechanics the micro-canonical con- 
dition rests on the further assumption of equipartition 
on the constant energy shell in the phase space, which 
in turn can be justified by the dynamical chaos caused 
by the non-linear interactions between particles. In this 
respect, the classical derivation is in better shape. In 
contrast, the quantum systems have no phase space, and 
also they are fundamentally linear. A typical state of 
an isolated quantum system is not an eigenstate but a 
superposition of eigenstates: 

N 

vf^^Cc/.,, (1) 

i=l 

where ^' is the wave function of the superposition, (pi is 
the wave function of the i-th eigenstate, Ci the corre- 
sponding complex amplitude, and N the total number of 
eigenstates. Therefore, the straightforward counterpart 
of the classical microcanonical condition would be to con- 
strain the possible choices of ^ to the "energy shell" in 



the Hilbert space: 

N 

Y^E^p^^E,,, (2) 

i=l 

where pi = |Cip are the occupation numbers of quantum 
states and i?av is the energy of the quantum superposition 
set externally and referred to below as "average energy." 

Condition ^ is, however, different from the conven- 
tional microcanonical condition, because the latter in- 
volves the summation only over the eigenstates inside a 
very small energy window ii'av ± 5E. 

Why the system should limit itself to a small energy 
window is difficult to justify unless, for example, one as- 
sumes that the quantum system is subjected to an exter- 
nal source of decoherence with the subsequent collapse 
of the density matrix. However, the introduction of col- 
lapse would imply that non-relativistic quantum mechan- 
ics is not self-contained, when it comes to justifying the 
Boltzmann-Gibbs equilibrium. 

Besides the conceptual issues, there are also practical 
ones. When a well isolated quantum system having not 
too many particles but many quantum levels is shaken 
in an experiment and then left to itself, the energy win- 
dow of the participating eigenstates can easily become 
larger than the temperature. Would such a system end 
up exhibiting Boltzmann-Gibbs statistics? 

It is clear a priori, that a significant departure from the 
narrow-energy- window constraint can easily lead to devi- 
ations from the Boltzmann-Gibbs statistics One can, 
consider, for example the case of two narrow energy win- 
dows. Still one can hope that somehow the "most proba- 
ble departure" from the narrow energy window condition 
would still support the Boltzmann-Gibbs equilibrium. 
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If one is to begin addressing the above issues, the un- 
avoidable hmit to consider is the system of >> 1 quan- 
tum levels with constraint ([2]) and no limit on the en- 
ergy window, i.e. all quantum superpositions of form ([T]) 
satisfying condition ([2]) are equally probable and, there- 
fore, the probability density as a function of the complex 
amplitudes Ci is proportional to the volume element on 
a manifold in the Hilbert space constrained by Eq.(I2]). 
Following Ref.Q, I call this condition "quantum micro- 
canonical" (QMC) to contrast it with the conventional 
micro-canonical condition. 

The general approach of assigning the probability on 
the basis of volume in Hilbert space has received a good 
degree of attention in recent years. Some of the rele- 
vant works [1, U, [E, Q applied this approach to the con- 
ventional micro-canonical case with small energy window 
for the participating eigenstates. Other works [3, 0, H, [1] , 
however, have looked precisely at the QMC alternative. 

In particular, it was found in Ref. Q for the case of 
equally spaced levels (and confirmed in the present work 
for the general case) that, as iV — oo, the volume of 
the Hilbert space as function of i?av acquires the charac- 
ter of a (5-function with peak located at E^v = X^i ^i- 
This result, however, does not imply that it is not impor- 
tant to consider the case of i?av different from the above 
value. The situation here is analogous to the conventional 
micro-canonical description, when the most probable po- 
sition of the narrow energy window would correspond to 
the infinite temperature, but one would still like to know 
the result for a finite temperature. 

The goal of the present work is to obtain from the 
QMC condition the statistics for the occupation num- 
bers of individual quantum states both for the entire 
isolated quantum system (Section |TT|, and for the den- 
sity matrix of a small part of it (Section IIIip . It is to 
be shown analytically that this statistics is dramatically 
different from the Boltzmann-Gibbs statistics: the oc- 
cupation numbers of quantum states decay with energy 
algebraically rather than exponentially, and, in addition, 
the macroscopic limit is routinely accompanied by con- 
densation into the lowest energy state. 

Section IIVI contains preliminary numerical tests of 
some of the above analytical results for finite systems. 
This section also reports a numerical finding that the 
expansion of the basis states of large but finite random 
matrix Hamiltonians in terms of the eigenstates of these 
Hamiltonians obeys the statistics found in this work. 

The implications of the above results are to be dis- 
cussed in the concluding remarksf Sect ion |V|) . 

The notion of chaos does not play any role in the forth- 
coming derivation, but it will also be touched briefly in 
the concluding remarks. 

Even though the formal treatment below is based es- 
sentially on the geometrical analysis of many-dimensional 
manifolds, it should not escape the readers that the 
end result is similar to the Boltzmann-Gibbs distribu- 
tion for grand canonical ensemble — though not for the 
average values of the occupation numbers (pi), which 



the Boltzmann-Gibbs distribution aims at predicting, 
but rather for the probabilities of variable pi to admit 
different values. In other words, it is an example of 
superstatistics[l0] — consequence of the fact that the oc- 
cupation numbers pi, which are interpreted as quantum 
probabilities, are themselves subject to the probability 
distribution. 



II. STATISTICS FOR THE EIGENSTATES OF 
AN ISOLATED QUANTUM SYSTEM 

A. Formulation of the problem 

The Hilbert space is parameterizable by the absolute 
values \Ci\ and the phases tpi of the complex amplitudes 
Ci = \Ci\e^'^'. Without the energy and the normalization 
constraints, the volume element in this space is given by 
dV = nf ICI d^, d\a\ = iJlf rf<^. rf(iaP). Both en- 
ergy and the normalization constraints involve only |Cip, 
and, therefore, when assigning the probabilities in the 
Hilbert space, one can integrate over phases ipi and then 
deal only with the subspace of occupation numbers 

P^ = |C,|2. (3) 

The Hilbert space volume element is then given by 

N 

dV^Y[dp, (4) 

i 

with an unimportant prefactor. It is always to be as- 
sumed that >> 1. 

In the rest of this paper, I will be calculating the vol- 
umes of manifolds in the Euclidean space of variables pi 
constrained by 

(i) energy condition ([5]); 

(ii) normalization condition 

N 
i 

and (iii) positivity condition 

P^ > 0, Vi. (6) 

Energies Ei are ordered by their values with the min- 
imum one being -Emin = Ei and the maximum one 
i^niax = En. The zero reference point for the energies 
is chosen such that 

N 

Y^E.^O. (7) 

i 

When not stated otherwise, it will be assumed below that 

Eav < 0. (8) 
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For a given value of index k, I denote the {N — 1)- 
dimensional Euclidean space of all variables {pi} with 
i^k 'AS {pi}k- 

The probability of the fc-th occupation number to have 
certain value pk is proportional to the volume Vk{pk) of 
the (N-3)-dimensional manifold in the {pij^-space — to 
be denoted as Mk — constrained by conditions (O in 
combination with 



N 



^ = 1 -pk, 



(9) 



-consequence of ([S]), and 



N 



(10) 



Condition pop . while obviously originating from re- 
quires a preliminary manipulation equivalent to shifting 
the origin of the energy axis to -Eav Namely, E^v in 
the right-hand side of Eq.([2]) has to be multiplied by 
Pi [equal to 1 according to (O] and then the result 
transformed to Ea. (|10p . Important for the subsequent 
derivation is the fact that energy hyperplane represented 
by Eq. pUj) crosses the origin of the {pi} ^-space, when 
Pk=0. 

The manifold Mk has a character of {N — 3)- 
dimcnsional polygon with flat faces, edges, etc., because 
all conditions constraining it represent hyperplanes in the 
{Pi}fe-space. 

The probability distribution of pk is then 
P{Pk) = Vkipk)/ (^Jo Vk{p'i,)dp',^^, and the average 
value of Pk is 



(Pk) = 



lo P'kVk{p'k)dp'k 
IoVM)dp'k 



(11) 



B. 



The case of pfc <C 1 



I first consider the case pfc ^ 1. 

The manifold constrained by conditions ^ [51 [T(I| can 
now be described as follows: 

The intersection of ^ and Q is a many-dimensional 
analog of a tetrahedron. It has dimension (iV— 2) — equal 
to that of the normalization hyperplane (O with (N — 1) 
vertices located at the intersections of the {N — 1) axes 
of the {pi}k space with the hyperplane i.e. in the 
{Pi}k-s>pa.ce, each of the vertices has coordinates of type 
(0, 0, 1 — Pk, 0) — all projections are zero, except 
for one, which is equal to 1 — pk- I call the resulting 
object "Hypertetrahedron" . This Hypertetrahedron is 
then cross-sected by the energy hyperplane (jlOp. 

When Pk is small, the renormalization of the volume 
Vk{pk) with respect to Vfe(O) can be decomposed into the 
"normalization factor" _FV due to the non-zero value of 



Pk in Eq. ([9]) and the "energy factor" Fe due to the non- 
zero value of Pk in Eq. pU]) : 



VkiPk)=VkiO)FNFE 



(12) 



The normalization factor is given exactly by 

FN^{l~Pk)''^^ (13) 

for large or small pk ■ It is the consequence of the fact that 
the change of {1 — pk) in the right-hand side of Eq.® 
rescales the distance between any point of the Hyper- 
tetrahedron and the origin of the {pi}k space by factor 
(l—pk). Since the energy hyperplane PH)) passes through 
the origin (at pk = 0), each of the {N — 3) dimensions of 
the intersection manifold simply undergoes rescaling by 
factor {1 — Pk) thus leading to factor (|13p . 

The calculation of the energy factor Fe requires more 
effort. The change of ~{Ek — Es,v)Pk in the right-hand 
side of Eq. (flO)) shifts the energy hyperplane in the trans- 
verse direction, but the resulting change of manifold Mk 
does not any longer amount to a self-similar rescaling. 

The volume 14 of the {N — 3)-dimensional manifold 
Mk can in general be presented as a product of {N — 3) 
characteristic linear parameters rjka '■ 



N-3 



(14) 



These parameters can be defined itcratively in the fol- 
lowing way: riki = Vk/Vk,N-i, where Vk,N-i is the 
volume of one of the {N — 4)-dimensional faces of Mk] 
Vk2 = Vk,N-4:/Vk,N-5, where where Vfc,jv-5 is the volume 
of one of the {N — 5)-dimensional faces of the {N — 4)- 
dimensional face selected in the previous step; etc. 

After the small shift of the energy hyperplane (fTU)) 
by —{Ek — Eiiv)pk, each linear parameter rjka changes 
slightly to 

VkaiPk) = r]ka{0)[l ~ \ka{Ek - £^av)Pfe], (15) 

where Xka are unknown rescaling coefficients. As a result, 

N--3 



Fe= l[[l-Xko.{Ek~E,,)pk] 



where 



-{N-3)\kiEk-Ei,-v)Pk 



(16) 



N-3 



N 



(17) 



Coefficients Xka are not well differentiable with respect 
to Pk and i?av, because the change of pk and i?av is ac- 
companied by the change in the number of vertices of 
manifold Mk- However, the internal self-consistency of 
the present treatment indicates, that the overall renor- 
malization factor Fe depends on pk and i?av sufficiently 
weakly and can be efficiently approximated. 
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Even though each renormahzation factor [1 — \ka {Ek — 
Ea,v)Pk] in Eq.((Tn]is very close to 1, the product of the 
(iV — 3) of these factors may be significantly smaller than 
1 without compromising the validity of small-pfc approx- 
imation (flB]) for Fe alone. Yet, when pk <C 1, but both 
Fn and Fe depart significantly from one, one can worry 
that the effects of shifting the normalization and the en- 
ergy hyperplanes and (jlOp ] do not commute with 
each other, and therefore, the resulting renormahzation 
is not equal to the product of F/v and Fe- This is, how- 
ever, not the case, because the shift of the normalization 
hyperplane amounts to a simple rescaling, and after that, 
the shift of the energy hyperplane always begins from the 
manifold of the same geometry. 

Central to the present work is the result that in the 
leading order in is simply independent of k. I de- 

note this independent value as A without a subscript. It is 
shown in Appendix \X\ that the linear parameters in (jl4p 
can always be chosen such that all but one summands 
are equal to each other in the expressions for two dif- 
ferent renormahzation coefhcients = ^ 

and A; = Y.a ^^ ^la' ^ka = for aU a except 
for one value ao- In a typical case, however, Xkao s-^d 
Xiao much smaller than the rest of their respective 
sums. 

Substituting A instead of Afe in (|16p and then combin- 
ing in (jl2p the resulting expression for Fe with Fn from 
while keeping only the leading order in N, I obtain 



VkiPk) = V^fc(0)e-^f^[i+^(^''-^-)], (18) 



As long as 



i + A(£;fc-£;av)» -, (19) 



Vk{pk) decays almost completely, when ^ 1, and, 
therefore, expression (fT8|) is sufficient to calculate (pk) 
from Eq. pTjl . which gives 



N[l + X{Ek-E,,)y 



(20) 



The value of A can now be found numerically from either 
of the following two conditions originating, respectively, 
from Eqs.® and dH): 



N 



k=l 



N 



Y,{Ek-E,,){pk) = 0. 



(21) 



(22) 



fe=i 



Expression ([^0]) for {pk) has the property that, if the 
value of A is found from one of the two conditions — 
(PTjl or ([m, the other one is fulfilled automatically. 

The value of A thus obtained becomes a function of 
average energy A[i?av]- Below I use A both with and 



without its argument. In order to distinguish the argu- 
ment of function A[-Bav] from the multiplication of A by 
an expression in parentheses, the argument of A[i?av]j if 
present, will always follow A in square brackets. 

It is useful to present the conditions (^11 also in 
the integral form with the values of (p^) substituted from 



1 

TV 



+ 00 



v{E)dE 



l + X{E- 



= 1; 



{E - E,^)v{E)dE 
\ + X{E- £;av) 



(23) 



(24) 



where {v{E) is the density of states corresponding to 
the energy spectrum {Ek} and satisfying the condition 
C^.{E)^N. 



C. Meaning of A 

The parameter A or, more precisely, A^A has the mean- 
ing of inverse Hilbert space temperature. It was intro- 
duced to describe the volume change of manifold Mk in 
the (A — l)-dimensional {pi}fe-space in response to the 
change in the right-hand side of the energy constraint, 
but, in the leading order in 1/N, it also describes the 
change of volume T4ot of the entire energy manifold con- 
strained by conditions (|2I5I6P in the full A^-dimensional 
Hilbert space of the problem as a function of Ea^: 

A[i?av] = ^a|;^l0g^tot(i?av) + 0(1/A^) (25) 

[The extra dimension would introduce only one extra lin- 
ear parameter 770 and one more coefficient Ao,q. in the sum 
of (A^— 1) other comparable coefficients in the early proof 
that justified the single value of A for all {p^jfe-spaces — 
see Section Hi Bl and Appendix [XI] As a consequence, 



Vtot(£^av) = KiaxBXp 



A^ 



X[E]dE 



(26) 



where Knax is the maximum value of Vtot corresponding 
to -Eav = 0, which, in turn, is the average value of all 
energies in the spectrum as defined by Eq.Q. 

In order to prove that the maximum of Vtot is indeed 
located at i?av = 0, one should note that, according to 
Eq. (p5)) . this maximum implies A = 0. Equation (|20|) 
then gives {pk) = 1/A, which, according to Eqs. (|7l22p . 
can only be the case, when i?av — 0. This general result 
is in agreement with the analysis of Refs.[0]. 

In principle, the point i?av = may or may not coincide 
with the maximum of v{E), which is already a significant 
departure from the conventional statistics predicting the 
most probable state of the system (zero inverse temper- 
ature) always at the maximum of i^iE). 
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FIG. 1: (Color online) Sketch of X[E]. Here E is the average 
energy identical with E^v in the right-hand-side of Eq.([2ll. 
Critical average energy Ec is defined by Eq. (f30)l . 

Another important difference is that even when the 
maxima of Vtot(-B) and i^iE) coincide, Vtot{E) decays ex- 
ponentially faster than v{E), which has the consequence 
that the small-pfe condition ((T9)) can be easily violated 
for the low-lying levels leading to a sort of condensation. 

A typical dependence of parameter A on the value of 
the average energy is sketched in Fig. [T] One should, in 
particular, remember that i?av < (the default assump- 
tion for most of this paper) corresponds to A > and 
vice versa. 



D. Beyond small 

In general, whether or not the preceeding description 
is sufficient for calculating the occupation of all quantum 
levels depends on the spectrum of the problem and on 
the value of i?av One should attempt to solve the above 
system of equations and see whether {p^) <C 1 for all k. 
If not, then one should use the following results extended 
to the case of {pk) ^ 1- 

I 



14(pfc) = ^fe(0)exp<^ (iV-3) 



where X[E] is to be determined self-consistently by solv- 
ing Eq. (|2T|) or ((22)) , though this time not just for a single 
value of i?av but rather for i?av spanning most of the al- 
lowed interval [i?min, ^-max] as required by the integral in 
([3T|) . A possible algorithm of the overall self-consistent 
solution in discussed in Appendix [Bl I chose not to ap- 
proximate (TV — 3) with TV in Eq.([3T|). because it appears 
to be an important correction for finite- A'' systems. 



Expression ([20)) has a pole as a function of E^, which 
I denote as E^. 

Ex^E^,-'^. (27) 

Parameter E\ in the present treatment is analogous to 
the chemical potential in conventional thermodynamics. 
Like A[ii^av], it is the function of the average energy, 
Ex[Ea,-v], and, likewise, I will be using square brackets 
to refer to the argument of this function. 

Condition (fTO)) is satisfied for all levels, when 

g.nin-i?A» ^^^^^\ (28) 

which for all practical purposes translates into 

Ex < S„un - Oil/N) (29) 
When A is positive and small, it corresponds to negative 
Eav sufhciently close to zero, and therefore, the preceed- 
ing solution is valid. However, it always becomes violated 
as soon as i5av departs significantly from zero. Below I 
will be using variable Ec to refer to the critical value of 
i?av corresponding to: 

Ex[Ec] = E,^in- (30) 

(see Fig. [T]) 

In order to understand the regime of large pk, it is 
necessary to appreciate, that once condition (fTO)) or (|^5)) 
is violated, A would continue to describe the response of 
Vk{pk) to the small shift of the energy hyperplane (fTO)) 
around pk = 0- However, the whole function Vk{pk) 
given by Ea. (fT5)) either decays too slowly or increases, 
and, therefore, the linear approximation for the power in 
(fT8)) becomes insufficient and Ea. (|20)) is not justified any 
longer. 

It is, however, shown in Appendix |B] that the deriva- 
tion of (fT8)) is amenable to the case of arbitrary p^. The 
result is 



(31) 

I 

It is possible to anticipate the outcome of the above 
self-consistent solution qualitatively. Once the condition 
(1^5)) is violated, the average occupation number of the 
levels affected becomes a significant fraction of one. The 
overall normalization constraint (|2ip then implies that 
the number of these exceptional levels should be small, 
certainly much smaller than N. Therefore, as £'av contin- 
ues decreasing beyond Ec, Ex may increase above E'min 



log(l-pfc)+ / X[E]dE 
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but will stay very close to i^min- This leads to an impor- 
tant approximation, namely: for i?av < -E'cj 



1 



(32) 



which one obtains by substituting E\ 
Eq.dail). 

The value of \Ec\ is comparable to |ii^min|, when N is 
large, but not exponentially large, as might be the case 
in the numerical studies (see below) and also in nano- 
sized systems having not too large number of particles, 
but very large number of levels. In these cases, the full 
calculation of the resulting statistics has to be done nu- 
merically. 

However, for a physical system, having macro- 
scopic number of weakly interacting components, further 
progress can be made analytically. 



E. Macroscopic system with non-degenerate 
ground state 

1. Definitions and assumptions about the macroscopic 
system 

As macroscopic, I understand a system consisting of 
macroscopic number Ns 10^^ of relatively weakly in- 
teracting parts. In a gas, one molecule would constitute 
such one part. In condensed matter systems with finite- 
range interactions, a part would imply a cluster of atoms, 
whose volume energy is much greater than the surface 
energy. Each part is assumed to be characterized by a 
finite Hilbert space. (The limit of large Hilbert space per 
constituent part is considered in the next subsection.) 
The total number of levels in such a macroscopic system 
satisfies inequality: 



iV > 2^= > Ns 



(33) 



Due to the large number of weakly interacting parts, the 
density of states I'iE) of this macroscopic system is as- 
sumed to have narrow Gaussian peak around E = 0, 
which is set by Eq.(I7]) to be equal to the average of all 
energies {Ei}. The mean-squared deviation of ^{E) from 
the above peak position is 



1 ^ 
as = -yE 



Nsao, 



(34) 



i=l 



where ctq is the typical mean-squared deviation for a con- 
stituent part. 

The average energy for each constituent part can also 
be set at zero. With the above convention. 



(35) 



where, -Emino is the typical minimum energy for a con- 
stituent part. 



In the conventional micro-canonical formulation of sta- 
tistical physics, the temperature T corresponding to en- 
ergy Eav is defined (with the Boltzmann constant set to 
l)by: 



1 _ d log iy{E) 
T ~ 



dE 



(36) 



E—Es, 



The substitution of Gaussian approximation for iy{E) 
then gives T — —OsjEa^,. (Positive temperatures cor- 
respond to £^av < 0.) In a typical situation of physical 
interest, T Eq^ where Eq is a characteristic one-particle 
energy in the Hamiltonian. Therefore, 



(37) 



Like in Eq. (j35p . the above inequality is the consequence 
of Ns ^ VN^- Therefore, it obviously extends to all re- 
alistic cases of relatively large constituent parts at tem- 
peratures in the range 10~^ — 10^ times Eq. 

It is shown in Appendix [C] and further in Appendix [D] 
that, under the above conditions, the critical average en- 
ergy Eq defined by Eq. (|50)) satisfies inequality 



\Ec\ < CTs < \En 



\E,: 



(38) 



Finally, important for proving the condensation into a 
single lowest level is the property 



Eo — E„ 



> 



\E^ 



1 



N 



NX[E, 



(39) 



where E2 is the energy of the second lowest level. This 
inequality is the consequence of the exponential small- 
ness of Ns in comparison to N in combination with the 
fact that E2 — i?min is, crudely speaking, a single particle 
property falling on the scale of |i?min|/A^s multiplied, per- 
haps, by some other factors depending polynomially on 
Ns- The rightmost expression in ([55]) is the consequence 
of ([5^ . Inequahty ([55)1 is illustrated in Appendix [Cl 



2. Results 

It is shown in Appendix ID] that . when condition (j39p is 
satisfied, E\ is pinned between i?niin and E2, sufficiently 
far from E2, so that 



A^[l + \{E2 ~ Sav)] = N 



E2 — Ex 

E\ — Eas 



> 1. 



(40) 



This justifies the approximation pk ^ 1 for k > 2 and, 
therefore, the validity of formula (pOI for the second low- 
est level and all levels above it. In this formula the 
value of A can then be very accurately approximated by 
Eq.®. 

As far as the volume Vi(pi) is concerned, it is narrowly 
peaked around the maximum, which simultaneously be- 
comes the average value of pi : 



(Pi) 



E, 



En 



(41) 
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Such a condensation into the lowest energy state amounts 
to a significant departure from the result of the conven- 
tional microcanonical recipe. 

In retrospect, it is also clear that for the case, when 
the occupation of only one lowest-energy level violates the 
condition (pi) <C 1, and therefore formula does not 
describe Vi{pi), still formula ([^0]) would give an excellent 
approximation for all (pk) including (pi), if the value of A 
in that formula is found self-consistently from Eq. ([2T|) or 
([22|) . The reason is that since formula (f20|) is supposed 
to describe accurately the occupations of all levels begin- 
ning from the second, the occupation of the remaining 
(first) level is bound by normalization constraint (|2ip to 
have the right value. The self-consistent solution using 
formula (PO)) would produce E\ coming as close from be- 
low to ii^min as ncccssary in order to reproduce the value 
(HT|) . As far as other levels are concerned, for the ab- 
solute majority of them, the approximation E\ « -Emin 
would remain very accurate independently of whether Ex 
is slightly above or slightly below -Emin- [Here, a few low- 
lying levels may constitute a possible exception related 
to the fact that there is some uncertainty in the present 
derivation about whether Ex stays much closer to £^min 
than to i?2- If it does, which I suppose is the case, then 
the above procedure would be very accurate for all levels 
with k > 2.] 



F. The limit of large number of quantum states 
per particle 

When a macroscopic system consists of particles hav- 
ing translational degrees of freedom, the kinetic energy of 
the particles can reach very high values, before the par- 
ticles are able to escape from the system. Therefore, the 
number of quantum states per particle in such a system 
can be very large. 

In order to analyse this limit in the simplest case, one 
can consider a system of Ng identical non-interacting os- 
cillators having energy levels equally spaced by ft and 
the average energy per oscillator ne^i, where rig is a fi- 
nite number. 

Each oscillator can, in turn, be described as a large 
spin S in magnetic field in the limit — *■ oo. The energy 
of this spin would he E — flSz, where the projection 5*2 
admits (25' -I- 1) values between —S. 

In this case, the ground state energy of the whole sys- 
tem is 



while the average energy is 

Therefore, according to formula (|¥T|). 

S -Ue, 



{pi 



S s- 



1; 



(42) 



(43) 



(44) 



In other words, the most probable state of such a system 
is the ground state with vanishingly small corrections — 
quite a surprising result. 

In order to understand it intuitively, one needs to re- 
member that the exact value of (pi) remains less than one 
[see Eq. (|D16|) ]. and, moreover, 1 — (pi) ^ If any 

eigenstate of the spectrum remains completely unoccu- 
pied on average, it means that the corresponding volume 
in the Hilbert space is zero. Therefore, each of many 
eigenstates above i?av has to have some non-zero average 
occupation. At the same time, the eigenstates below iJav 
need to have much greater occupation in order to balance 
in Eq. ([22|) many more eigenstates above ii'av It simply 
turns out that the volume of the Hilbert space is max- 
imized, when almost all (but not all) of the probability 
weight goes into the ground state. 



G. Typical pure state 

Even though the statistics derived so far has been ob- 
tained through averaging over all possible quantum states 
subject to the QMC condition, the resulting statistics 
also describes a typical one among them in the following 
sense. Once a single state is selected, it will have very 
large number of eigenstates in each small energy interval 
between i^niin and -Emax- Individually, the occupation 
numbers of these eigenstates will fluctuate according to 
the probability distribution proportional to their respec- 
tive Vk{pk)- That distribution will depend only on the 
energy of each of these eigenstates, and therefore, within 
a small energy interval, it will be approximately the same 
for all of them. Consequently, the average occupation 
number of eigenstates within any small energy interval 
will be given by formulas for {p^) obtained above. 



III. SMALL SUBSYSTEM WITHIN A LARGE 
ISOLATED SYSTEM 

A. Formulation of the problem 

Now I proceed with deriving the energy distribution 
for a subsystem of an isolated system — subject to the 
QMC condition. It is assumed that the subsystem and 
the rest of the system — environment — do not interact 
with each other. Therefore, the eigenstates of the whole 
isolated system can now be labelled by two indices as 
follows: 



(45) 



where indices a and (3 and the corresponding eigenstates 
ipa and (jjp refer to the subsystem and the environment 
respectively. The subsystem has Ni states with energies 
Esa ■ The environment has N2 states {N2 ^1) with 
energies EsfS- The zero reference point for each set of 
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energies is chosen such that 



Esc = 0, 



and 



N2 



= 0. 



(46) 



(47) 



The energy of each eigenstate '^ap of the whole system 
is then 



Eap — Esa + EeP- 



(48) 



The occupation number of each eigenstate is Pafi- In 
this formulation, the density matrix of the subsystem, 
denoted as psaa' , has only diagonal elements. 
I now focus on finding the diagonal element 



N2 

PSaa = Pa — ^ Paf3- 
13=1 



(49) 



Here I defined variable pa just to shorten the notation. 



The total normalization constraint and the energy con- 
straints now have form, respectively: 



and 



^Pa+^Pb = 1, 



EaPa + ^ EbPb = -Ea 



(52) 



(53) 



Given Eg. ([5^ . condition ([53]) can be replaced with 

Y,iEa ~ E,,)pa + Y^^Eh - E,^)pb = 0. (54) 



The a- and the 6-states can now divide between them- 
selves the occupations and the total energy as follows: 



Y(.Ea - E,^)pa = Ea, 



(55) 



(56) 



B. General solution 

I now re-label the N2 states contributing to pa with 
index a, and the remaining (A^i — 1)N2 states with index 
b. This results in two new sets of occupation numbers 
and energies: (j)a,Ea) and (pf,,i?f,). Subscripts a and b 
will play dual role below: as indices and as labels of two 
different sets. The summation over a implies the first 
set, and a summation over b implies the second set. In 
cases, when I refer to the individual members of each set, 
I use the "label-and- number" subscript such as, e.g., Ea2, 
which refers to the second lowest energy of the a-set, or 
Ebmin refers to the minimum energy of the 6-set. 

The new sets of energies have the average values, re- 
spectively: 



^ N-2 

a— 1 



N2{Ni - 1) 



E. 



Sa 



b=l 



Ni - 1 



(50) 



(51) 



YPb = 1 - Pa> 



Y,{Eb - E,^)pb - -Ea, 



(57) 



(58) 



where Ea is the difference between the average energy 
of the a-set and i?av It is an auxiliary parameter to be 
determined simultaneously with pa ■ 

The goal now is to obtain {pa) — the average value 
of Pa over all points in the Hilbert space constrained by 
conditions (l52l [54l) in combination with 



Pa,Pb > 0. 



(59) 



The probability of each pair of values {pa, Ea) is pro- 
portional to the volume in the Hilbert space constrained 
by conditions (|55ll59p . The constraints on a- and 6-states 
can then be treated independently by analogy with the 
problem for the whole system that led to Eg. (I3T|) . This 
results in the following expression: 



V{pa,EA) = Voexp I N2\ogPa + ^2(iVi - l)log(l - Pa) + N2 



Xa[E]dE + N2{Ni - 1) 



1-Pq 



Xb[E]dEj 
(60) 
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where \a{E) and \b{E) are the parameters analogous to 
A introduced below for the spectrum of the entire sys- 
tem, but this time defined for the spectra of Ea and Eb 
separately, and Vb is an unimportant prefactor. 

Due to the fact that N2 ^ 1, the expression ([TO]) must 
be very sharply peaked near the maximum of the function 
in the power of the exponent. Therefore finding (p^) is 
reduced to finding the value of pa at the maximum of this 
power. In order to locate that maximum, I look for the 
zeros of the partial derivatives of the exponential power 
in (|60p with respect to Ea and pa- Differention with 
respect to Ea gives 



A^i - 1 

I - Pa 



E, 



Ea 



Pa 



Xa 



(61) 



Now, differentiating the power in Eq. (|60p with respect to 
Pa and also using ([6T|) . I obtain 



jOq(1 - Nipa) = Ea Xa 



E, 



Ea 

Pa 



(62) 



One can get a useful insight into the solution of 
Eqs. (l61|62|) by substituting Aa,Af, with new variables Ex^ , 
Ex,: 



Xa[E] = 



Xb[E] = 



1 



E-ExAE] 



1 

E-ExAE] 



(63) 



(64) 



Subscripts a and b in variables Aa,Af,, Ex^, Ex, are just 
the labels of the characteristics of the respective energy 
spectra, i.e. they are not indices running over a set of 
integer values. 

After some manipulations, the above substitution gen- 
erates two equations equivalent to (|61l62p : 



Ex 



E,, 



Pa 

Ea 

Pa 



E,. 



= Ex 



Ex 



E, 



Ea 

Pa 



Em — 



Ea 



1 



Pa 



(65) 



(66) 



The latter equation is the key to the following solution 
for the macroscopic environment. 

In general, Egs.lpTl [S^ . or equivalently Eqs. (|65l66p 
should be solved numerically. However, for the case of 
a small subsystem and a macroscopic environment, and 
with the realistic value of E^v for the whole system (as 
discussed in Section UTE]) . the system of Eas. (|65l66p can 
be solved analytically. 

There is also another analytically solvable limit, which 
corresponds to the case of high Hilbert space tempera- 
tures, i.e. very small Xa and Xb- This limit is not to be 
considered in this work. 



C. Subsystem in a macroscopic environment 

The condition of macroscopic environment and a small 
subsystem amounts formally to the presence of narrow 
Gaussian-like maximum around = in the density of 
states of the environment i>EiE) with the mean-squared 
spread of energies as satisfying the conditions: 

\EEyainl EEme,^, {E^l > y/^ > \Esa\, Va. (67) 

In addition, there is a reasonable condition for the differ- 
ences between two lowest energy states for the subsystem 
and the environment: 



£'52 — Esn 



> Ee2 — E 



i?min ■ 



(68) 



It is shown in Appendix [El that . in this case, the occu- 
pations of the lowest subsystem state in the leading order 

of Esa/EEmin is 



E,: 



EEn 



1 



E, 



Eeu 



(69) 



and, for the remaining states with a > 2, 



{pa) = 



1 



1 - 



E,, 



EEn 



(70) 



where the right-hand side is obviously independent of a. 

It is also possible to obtain more general formulas (see 
Appendix IE|): 



(Pi 



and, for a > 2, 



Err 



Ni-1 



(Pa) 



1 -Bav — E 



Emin 



Ni Esa - EEn 



(71) 



(72) 



where E'min — Esmin + EEmin- Formulas (|71l72p are cer- 
tainly valid up to the first order in Ega/ EEmim but in 
fact have a broader range of applicability, because con- 
ditions (|67l68p are sufficient but not necessary for the 
validity of approximation (j71l72p . For example, this ap- 
proximation also describes the case Esa ~ Esmin, when 
condition I^sqI *C l^'Bminl is replaced by the require- 
ment that TVi ^ 1 and the density of states for energies 
Esa have a Gaussian-like narrowly peaked shape around 
E — 0. Further discussion of the necessary conditions 
for the validity of approximation (|71l72p is given in Ap- 
pendix |e1 

The basic assumptions leading to the above results 
contain a loophole of neglecting the interaction between 
the subsystem and the environment. The same loop- 
hole is also present in the conventional micro-canonical 
derivation of the Boltzmann-Gibbs statistics. In prin- 
ciple, given the condensation of the entire system into 
the lowest energy state, one should not be surprised that 
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a similar property is exhibited by a subsystem. Yet 
one can still worry about the validity of the condensa- 
tion into the single lowest energy state of the subsystem 
[Eq. (j69p ]. when the interaction energy with the environ- 
ment is much greater than the separation between the 
lowest and the second lowest energy levels of the subsys- 
tem. In this case, the occupation numbers of the sub- 
system would depend on entanglement properties with 
the environment in the ground state of the whole system. 
This loophole potentially opens the window for chaos and 
non-integrability to play a role in the resulting statistics. 
It also cannot be excluded that the result may then repro- 
duce the Boltzmann-Gibbs statistics in the energy range 
of the order of the subsystem-environment interaction en- 
ergy. This issue is not addressed further in the present 
work. 



D. Typicality for the pure states of the whole 
system for the density matrix of a subsystem 

The probability distribution of parameters describ- 
ing the density matrix of a small subsystem of the 
whole system exhibits exponentially narrow l/A'^2) 
maxima controlling the average values of these param- 
eters. Therefore, a random choice of QMC-constrained 
single quantum state would be, with probability 1 — 
0(l/iV2), exponentially close to the average values com- 
puted above. This situation is analogous to the "canon- 
ical typicality" 0, ^ for the conventional microcanonical 
condition. 



IV. NUMERICAL TESTS 

The analytical results of this work and there relevance 
to realistic systems should be checked numerically. Here, 
I present two such tests of preliminary nature addressing 
only the statistics for the entire isolated quantum system. 



A. Direct random sampling 

The first test is a direct Monte-Carlo sampling of the 
Hilbert space under constraint ^ . I have done this sam- 
pling directly in the Eucledean space of variables {pi} us- 
ing the algorithm consisting of the following steps: (i) se- 
lection of an orthonormal basis in the (iV— 2)-dimensional 
hyperplane constrained by Eqs.(I21[5]); (ii) identification in 
that hyperplane a (A'^ — 2)-dimensional hypercube, which 
encloses all the vertices of the intersection manifold; (iii) 
random sampling of points within that hypercube; and, 
finally (iv) acceptance of only those random points, which 
in the original A'^-dimensional {pi}-space have all non- 
negative coordinates as required by constraint 

This algorithm is not very efhcient: its acceptance rate 
at step (iv) decreases by about factor of 10 as N increases 



(a) 

N=10; £av=-3. 
0.35 , — ' — ' — ' — , — ' — ' — ' — ' — , — ' — ' — 




0.00 ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — L 

-5 5 



(b) 

N=10; £av=-8. 




FIG. 2: (Color online) Two results of direct random sampling 
of quantum states in a system of ten energy levels. Average 
energies are indicated above the plots. Dots with error bars 
represent the averages over all sampled states. Solid lines 
represents the prediction of Eq. (|20p with the value of A ob- 
tained numerically from formula Ea. (l22|l . The inset in plot 
(b) magnifies the small {pk) part of the main plot. 

by one. Using Mathematica software, I was able to gener- 
ate statistically significant number of random points for 
the case of iV = 10 with the specrum and the average 
energy shown in Fig. [21 

This figure compares the average occupation numbers 
obtained numerically with the approximate theoretical 
values obtained on the basis of Eg. ipPl) with with the 
value of A found by solving Eq. (|22| numerically. 

The lowest level in Fig. [U^b) violates the condition 
(Pfc) ^ 1/-^ ^ 1- However, as discussed at the end of 
Section III El the overall structure of the more accurate 
solution guarantees that, for a single level violating the 
above condition, the result ((22l) would still amount to a 
very good approximation, even though the corresponding 
Hilbert space volume Vi (pi ) is does not any longer decay 
exponentially but instead is peaked around pi « (pi). 

Given that iV = 10 barely qualifies as a very large 
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number, the overall agreement exhibited in Fig. [2] is sur- 
prisingly good. 



B. Random matrices 

Now I show that the statistics obtained in this work 
manifests itself in a system described by a random matrix 
Hamiltonian. 

If one of the original non-eigenstates in the basis, where 
the random matrix is defined, has energy E^^ (the diago- 
nal element of the random matrix corresponding to that 
state), this imposes constraint ([2]) on the eigenstates par- 
ticipating in the expansion of the selected state. The 
present numerical experiment was based on a guess, suc- 
cessfully confirmed by the end result, that the eigenstates 
in this case would participate in the expansion of a non- 
eigenstate, as if that expansion was done randomly on 
the basis of the QMC condition. 

I took a 4096 x 4096 matrix, where all diagonal ele- 
ments and a fraction ^^|g of off-diagonal elements were 
assigned random values picked in the interval [—1,1]. 
The remaining off-diagonal elements were zeros. The 
Hamiltonian was diagonalized, and then one state of the 
original non-eigenbasis was chosen and expanded in the 
eigenstate basis. The weight of individual eigenstates 
fluctuated as expected from Eg. p^ . However, once the 
spectrum is divided in groups of 64 adjacent eigenstates 
having approximately the same energy, then the average 
weight within each group begins converging to the theo- 
retical approximation (PO)) as discussed in Section Hi Gl I 
further improve the error bars by combining the overall 
statistics for 34—35 non-eigenstates with average energies 
within a narrow energy window .Eav ± 0.01, where -Eav is 
equal to —0.5 and —0.9 in the two examples shown in 
Fig. [3l The good agreement is then revealed with the 
theoretical approximation (|20p . which uses A computed 
numerically from Eg. iP^ . 

When I increase the fraction of non-zero off-diagonal 
elements in the random matrix, the agreement between 
numerics and the theory continues to hold. In this case, 
however, the width of the nearly semi-circle eigenspec- 
trum increases, while the window for Eav determined by 
the diagonal elements of the Hamiltonian remains the 
same. As a result, A becomes small and the dependence 
(I20p becomes difficult to distinguish from a linear one. 
On the other hand, if the fraction of non-zero off-diagonal 
elements decreases, then the assumption of perfect mix- 
ing of eigenstates in each of the original basis states be- 
comes increasingly inadequate, as the weights of eigen- 
states start to peak around E = Eav Such a behavior 
is natural to expect as this system gradually approaches 
the limit of small off-diagonal elements in the Hamilto- 
nian (see, e.g., Ref.[ll|). 



(a) 



A'=4096; E^,=-0.5 
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W=4096; £'av=-0.9 
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FIG. 3: (Color online) Average occupations of eigenstates (p) 
participating in the expansions of non-eigenstates of the initial 
basis, which (the non-eigenstates) are selected from a narrow 
energy window iJav ± 0.01 around the values indicated above 
plots (a) and (b). Broken Unes in these two plots represent 
numerical results averaged over groups of 64 adjacent eigen- 
states. Smooth lines represents the prediction of Ea. (|20|l with 
the values of A obtained numerically from Ea. (|22|l . The den- 
sity of eigenstates z^(-E) for this system averaged over groups 
of 32 states is presented in plot (c) . 
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V. CONCLUDING REMARKS 

1) The results presented in this work, indicate that 
the statistical description of an isolated quantum system 
subject to a fixed energy constraint and unrestricted par- 
ticipation of eigenstates contradicts (at least in the limits 
considered) to the Boltzmann-Gibbs statistics derivable 
on the basis of the conventional microcanonical assump- 
tion. 

In particular, the resulting energy distributions for 
both the whole isolated system [Eq. ipD]) ] and a small sub- 
system of it [Eq. (f70|) ] show algebraic rather than expo- 
nential dependence on the energies of participating states 
as well as routine macroscopic occupancy (condensation) 
for the lowest-lying energy states — Eqs. (|41|44|69p . 

One should be mindful though of the loophole associ- 
ated with the neglected interaction between the subsys- 
tem and the environment — see the discussion at the end 
of Section Unci 

2) The statistics derived in this work is supported by 
the numerical findings presented in Section [IVI Particu- 
larly interesting is the finding presented in Section lIVBi 
that the expansion of the non-eigenstates of the random 
matrix basis in terms of eigenstates follows this statistics 
for large but finite random matrices, which are not too 
sparce. 

3) The finite size version this statistics might thus be 
observable after a strong (and preferably non-integrable) 
perturbation of a well isolated system having a relatively 
small number of particles but a relatively large number 
of quantum levels. Nanoscale-limited systems should be 
good candidates for such a study. For a better control of 
the total energy after the perturbation, one can proceed 
in analogy with the "numerical experiment" on a random 
matrix presented in Section [IV Bl Namely, one can force 
the system into a single quantum state, e.g. the ground 
state, before perturbing it. 

4) In systems of bosons, the condensation described 
in this work into the ground or a few lowest states may 
produce an appearance of Bose-Einstein condensation, 
because the ground state is indeed Bose-condensed. Yet 
the nature of the two kinds of condensations is different. 
The former represents a jump in the occupation of many- 
particle states, while the latter is a single-particle phe- 
nomenon accompanied by the usual exponential statistics 
for the occupations of all many-particle quantum states. 

5) When it comes to macroscopic systems, the contra- 
diction indicated in remark 1 reenforces the concern that 
non-relativistic quantum mechanics alone is not sufficient 
to justify the Boltzmann-Gibbs statistics. One needs an 
assumption of external origin, such as the quantum col- 
lapse of the broad distribution of eigenstates into a nar- 
row energy window postulated by the conventional mi- 
crocanonical description. 

If such a collapse happens even once, it appears very 
difficult if not impossible to realistically perturb a typical 
macroscopic system containing many weakly interacting 
parts into a state characterized by a broad energy range 



of participating eigenstates. 

Yet, if this collapse occurs continuously, its description 
would go beyond the linear quantum mechanics and in 
particular may imply additional source of energy fluctu- 
ations for the entire system. 

6) Even though the statistics obtained in this work 
appears to contradict to the everyday experience well 
describable by the Boltzmann-Gibbs statistics, it is still 
interesting to think, what the present statistics might im- 
ply, if one assumes that the entire Universe is describable 
by a single wave function. 

7) Quite a few researchers including this author (see 
e.g. Refs.[ll M, Q M, [H, [13, d) share the feeling 
that the elusive notion of quantum chaos plays an im- 
portant role in the foundations of quantum statistical 
physics. Chaos, however, plays no role in finding the 
most probable quantum state both in the present and 
in the conventional micro-canonical formulation. This 
suggests that the role of chaos is not to determine the 
equilibrium itself but rather to infiuence how a subsys- 
tem relaxes towards the equilibrium. 

Yet, one should not forget about the loophole related 
to the neglected interaction between the subsystem and 
environment [see the end of Section IIII C| . In addition, 
what quantum chaos certainly does is that it strongly 
suppresses the fluctuations of the number of energy lev- 
els within any fixed energy window — consequence of 
the repulsion of energy levels. Whether and how this 
property would affect the macroscopic characterization 
of equilibrium in quantum systems is not clear to this 
author. 

8) From a broader perspective, the statistics based on 
simple constraints ([21 [5l [6]) describes a distribution of an 
essentially positive and limited in the amount quantity p 
among N agents having characteristics {Ek}- It is, there- 
fore, tempting to speculate that such a statistics might 
be applicable beyond the quantum mechanical problems, 
in particular, to the problems of economics, when and if 
one finds a meaningful interpretation for constraint 
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APPENDIX A: PROOF THAT ALL VALUES OF 
Afc ARE APPROXIMATELY EQUAL 

In order to be specific, let us compare the intersection 
manifolds Mi and M2 pertaining respectively to Ai in 
the space of all variables pi excluding pi and A2 in the 
space of all variables pi excluding p2- Let us also, for this 
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part only, shift the origin of the energy zero point to i?av, 
i.e. the condition i?av = replaces Eq.([7]); and assume 
that El and E2 are not necessarily two lowest energies of 
the spectrum, but rather two arbitrary ones. 

In the first case, the intersection manifold is defined by 

P2,P3,-;PN >0, (Al) 
P2+P3 + ■■■+PN ^l, (A2) 

E2P2 + E-iP3 + ... + Enpn = V, (A3) 

where v — —Eipi. 

Parameter Ai characterizes the change of volume Vi of 
the above manifold in response to small shift v of hyper- 
plane (|A3|1 . Specifically, 

Ar=^^ =E-^=EAi.. (A4) 

^—^> Q — 1 ' a— 1 

In the second case, the intersection manifold is 

Pi,P3, ...,PAr > 0, (A5) 

Pi +P3 + ••• + = 1, (A6) 

Eivi + E-iPj, + ... + EnPN = V, (A7) 

where v — —E2P2- 

Parameter A2 characterizes the change of volume V2 of 
the above manifold: 

^^ = v-r =E-^ = Ea... (A8) 

If one changes the variable pi to p2 in Eas. (|A5IIA7[) . 
then one finds that the problems of calculating Ai and A2 
are nearly identical to each other with the only difference 
that term £'2^2 in Eq. (|A3p has to be replaced with Eip2 ■ 
One can already anticipate the result Ai ~ A2 from the 
fact that vectors {E2,E^, ...,Ejm} and {Ei.E^, ...,En} 
determining the normal directions of the hyperplanes 
(|A3P and (jA7|) . respectively, are nearly parallel to each 
other — consequence of the large number of identical 
components, the proof given below, however, takes a dif- 
ferent route. 

It is, in fact, possible to show, that the sets of charac- 
teristic linear parameters {jyia} and {7720} can be chosen 
such that rjia{v) — 7720 (w) for all but one value of a, 
and, therefore, — 4 out of A^ — 3 terms contributing 
to Al and A2 in Eqs. (|A41 lASp are identical to each other. 
Therefore, in general, 

A2 = Ai + 0(l/A^). (A9) 

The proof of the above statement lies in the fact that 
two {N — 3)-dimensional manifolds Mi and M2 (with 



the pi-to-p2 variable change in the second case) have a 
common {N — 4)-dimensional face. This face is defined 
by condition p2 = 0. Therefore, one can choose the first 
linear parameter t] in each case differently by dividing the 
volumes Vi and V2 by the (A^ — 4)-dimensional volume 
of that face, but then the rest of the parameters can be 
chosen identically, because they will describe the volume 
renormalization for the same face. 

I now describe the above (A^ — 4)-dimensional face in 
more detail and, in particular show, that the (A^ — 1)- 
dimensional hyperplane P2 = contains more than A — 3 
vertices of Mi or M2 — a necessary condition to form a 
{N — 4)-dimensional face. 

All the vertices of either manifold Mi or M2 are ob- 
tained by the intersection of the linear edges of the nor- 
malization Hypertetrahedron (defined in Section III B[) 
with the energy hyperplane. In turn, all of these edges 
have form pi — 1 — pj, (0 < Pi,Pj < 1) in the two- 
dimensional plane defined by condition = 0, Vfc 7^ z, j. 

The energy hyperplane (with v = 0) intersects such 
an edge only in the case, when Ei and Ej have opposite 
signs. (Here and everywhere below I ignore the non- 
generic case of Ei or Ej equal to zero.) Therefore, if 
there are K quantum states with Ei > Q [Ei > Eg,v] 
and L = iV - 1 - a: with E, < [E, < E^.-^], then the 
intersection manifold has KL vertices. 

Among these KL vertices, only those originating from 
the edges involving p2 change between the manifolds Mi 
and M2 (shift or disappear): 

if Ei,E2 > 0, then L vertices shift, possibly signifi- 
cantly; 

if i?i, -B2 < 0, then K vertices shift; 

if El and E2 have opposite signs, then L vertices in 
one case are replaced by K vertices in the other case. 

All other vertices remain identical. They all lie in the 
hyperplane P2 = 0. Their number is greater than KL — 
(A— 1), which, in turn is greater (normally, much greater) 
than (A^ — 3), when K,L > 2. Therefore this number is 
sufficient to form an (A — 4)-dimensional face. 

One can further show that the above vertices do not 
fall into a lower dimensional subspace but instead, can be 
used to form (A — 4) linearly independent vectors. The 
coordinates of these vectors can have (A^— 2) projections 
in the original Hilbert space: P3,P4, ■■■,Pn- Let us further 
assume that i?3 < and i?4 > 0. Therefore, there will 
be one vertex in the {^3,^4} two-dimensional plane. I 
denote this vertex as W34. In addition, there will be K' 
vertices in {p3,pi}-planes, such that i > 5 and Ei > 0; 
and L' vertices in {p4,pj}-planes, such that j > 5 and 
Ej < 0. The total number of vertices in the latter two 
sets is K'+L' — N—A. One can now form (A — 4) vectors 
by subtracting the coordinates of vertex w^^ from each of 
the above (A^ — 4) vertices. The resulting vectors will be 
linearly independent, because the projections on each of 
the axes p^^pQ, ...,pn will be present in this set only once 
in one of the vectors, and therefore the linear combination 
with other vectors cannot cancel that projection. 

Finally, the above {N — 4) dimensional face changes 
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identically under the shift of the energy hyperplanes in 
the both cases (which includes the possible change in the 
number of vertices). Indeed, because of the condition 
P2 = 0, the projected shifts of the energy hyperplanes in 
both cases are described by identical equation; 



£^3^3 + ... + EmPN = V. 



(AlO) 



Therefore, all ( A^— 4) renornialization coefficients describ- 
ing the volume change of this face with v can be chosen 
identically. 

As indicated in the discussion of manifold vertices, the 
structures of the two manifolds compared can be very 
different outside the common face. However, all these 
differences contribute to a single linear parameter 77 cor- 
responding to the direction perpendicular to the common 
face, which changes only weakly under the change of v. 

The last step for the complete rigor of the present 
proof would be to impose the limits on the exceptional 
cases, when the renormalization of a single parameter 
^kaa would be comparable with the sum of the rest of 

One possibility, of course, is that the latter sum ac- 
cidentally turns to zero, even though each term of that 
sum is comparable to Xkao ■ However, this would actually 
mean that the typical A, as used in the main text, is zero 
within the accuracy of the present approximation. 

The real problematic case would correspond to Xkao 
being much greater than almost all of other \ka- Such a 
situation appears to arise only in the case of little phys- 
ical interest, namely, when \Ek\ is of the order of the 
root-mean-squared value of all other energies, but this 
claim requires further proof. Here, I would only like to 
comment that the common face ( A^ — 4)-dimensional face 
is highly unlikely to be uncharacteristic of the (A^ — 3)- 
dimensional manifold, because, if all the vertices of this 
face are removed from the manifold, the remaining ones 
will not be enough to form an (A^ — 3)-dimensional man- 
ifold (when K > I and L> I). 



APPENDIX B: DERIVATION OF GENERAL 
FORMULA dUD FOR THE HILBERT SPACE 
VOLUME OF MANIFOLD Ah 



The derivation of formula (PT|) proceeds as follows. 
Both sides in each of Eqs.([ni[ini) are divided by {l~pk) 
and new variables 



Pi 



Pi 



l-Pk 



(Bl) 



for i k are introduced. This gives, respectively, 

N 



(B2) 



N 



where 



iji^k 



{Ek — Eg_^)pk 



l~Pk 



(B3) 



(B4) 



not a small number. The above equations in combination 
with constraints 



(B5) 



describe a manifold M'^. with volume VI. in the space of 
variables p ■ . Due to rescaling of axes (|Bip , 



Vk = {l-Pk 



(B6) 



where Vk is the volume of the original manifold Mk de- 
fined in Section III Al and power A^ — 3 is equal to the 
dimension of manifolds M'^ and Mk ■ 

Since the set of energy coefficients in the left-hand side 
of Eq. (jB3[) for manifold M(. is the same as the set for the 
original manifold Mk in Eq. pUj) , both manifolds are char- 
acterized by the identical dependence of the parameter A 
on the shift of the energy hyperplanes: 

1 d \ogV' 1 d \ogV' 
A[i?av] - ~ (B7) 

(see also Appendix [A|) . Therefore, 



(Af - 3) / \[E]dE \ 



(B8) 



where V^ is the volume of corresponding to w ~ 0, 
which is simultaneously equal to T4, when pk — 0. Com- 
bining Eqs.dm, (HSI and ([B4|, I obtain 



TT (-Bfc-Eav)Pfc 
Ea-v T-^, 



Vk{pk) = Vk{Q)ei^v\{N -i)[\og{l-Pk)+ \[E]dE\, (B9) 



I 

the same as Ea. (PT|) in the main text. This equation has to be solved together with Eq. ([^^ . which requires find- 
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ing X[E] for most of the allowed interval [i^min, -Bmax]- In 
practice, however, if the purpose is to find Vk{pk) for a 
particular value of i?av < 0, the interval [i?av,0] would 
surfice. This self-consistent solution can proceed as fol- 
lows: 

One first finds the values of A[-Bav] hi the interval, 
where the approximation (|20[) is valid, i.e. for ii'av be- 
tween and a certain value somewhat below the critical 
value Ec defined by Eg. ([50]) , and then proceed with 
reducing i^av further in sufficiently small steps and us- 
ing the approximation (|20p only for Ek > i?av, while for 
Ek < -Eav using the full formula (jB9p . where at each step 
the integral would require only the knowledge of X[E] 
for E > Ei^y. In principle, that integral may extend to 
-E > 0, but this corresponds to sufficiently large values of 
Pfc, for which Vk{pk) is guaranteed to exhibit fast expo- 
nential decay. Therefore, a cuttoff can be imposed 
that does not allow the upper integration limit to extend 
above 0. This cutoff is further discussed in Appendix ID] 



APPENDIX C: SYSTEM OF MANY SPINS 1/2 

In order to appreciate certain general aspects of the 
macroscopic case, it is sufficient to consider an otherwise 
very artificial example of Ns ^ 1 non-interacting spins 
1 /2 in magnetic field. The number of levels in this system 



(CI) 



In the basis of spins quantized along the direction of 
the magnetic field (z-direction) , each eigenstate is deter- 
mined by a set of spin projections Snz = ±1/2, and the 
corresponding energy is 



(C2) 



where Q, is the "Larmor energy" associated with the split- 
ting of spin states. The minimum and the maximum 
energies of this spectrum are respectively: 



and 



-Einax — ^Na^. 



(C3) 



(C4) 



The average of all energies in the spectrum is zero, and 
the mean-squared deviation from this average is 



All energies have form 



E = Enun + ^m, 



(C5) 



(C6) 



where m is an integer number between and Ns- The 
degeneracy of the energy levels corresponding to a 
given value of m is 

NJ 



m\(Ns — my. 

For m, (Ns — m) ^ 1, one can use the Stirling approxima- 
tion ml « -\/27rme™^^°S'"-i) and likewise for (TVs — w)! 
and Ngl, and then divide the degeneracy Dm by distance 
between levels to obtain the density of states: 



1 



^Afjog^^^+mlog^ 



n y 27rm{Ns - m) 

(C8) 

where m = [E — i?,iiin)/r2 — as follows from (|C6|) . The 
density of states can be further approximated by Gaus- 
sian: 



N 



(C9) 



which is the consequence of the Central Limit Theorem. 
It becomes increasingly inaccurate for |E^| 2> ^Jo^. 

As the most representative physical situation, let us 
consider the case T = fi, where T is the usual tem- 
perature defined by Eg. ([5^ . In this case, the Gaussian 
approximation (|C9|) for would imply Eav = —\^Ns = 
Emin/'2. This approximation is not very accurate for 
-Bav ~ Einin. For example, the accurate approximation 
Jell gives T = f]/log3 for Eav = Emi„/2. Yet this dis- 
crepancy does not invalidate the general estimate that in 
the physical range of interest 



lEavl > 



(CIO) 



with most routine being the case when i?av is negative 
and a finite fraction of i?min- 

Now, I would like to show that condition (|29p. which 
guarantees the smallness of all pfc, is violated for the 
present spin system almost immediately after E^v starts 
decreasing below zero. Let us consider what it takes to 
satisfy Eq.(lM]), when condition is fulfilled. The inte- 
gral in this equation has a character of the average value 
of {E — Eg^-y) under the effective distribution described by 
a product of slow varying function 1/ [1 + A(i? — i?av)] and 
sharply peaked symmetric function ^{E), which for the 
present purpose is well approximated by the Gaussian 
()C9p . In the leading order, this effective distribution will 
remain symmetric but with respect to a maximum, which 
is slightly shifted relatively to that of v{E). Eguation ([M|) 
can only be satisfied, when that maximum coincides with 
Escv. Such a condition gives A = E^^jos. Then requir- 
ing that E\ 

denoted earher [Eg. (f5D)l ] as Ec 



Einin, I obtain the critical value of i?av 



E, 



1 



c 



Eri 



-n < as, 



(Cll) 



which implies that the small-p^, condition (|29|) is violated 
at least for some fc, when the values of Eav fall in the 
range of primary physical interest (jClOp . 



16 



The analysis in Section III El relies on the assumption 
that, for iV, ~ 10^3, 



E2 — Eniin ^ 



\Er. 



N 



(C12) 



where E2 is the energy of the second lowest level, i?av sat- 
isfies (jClOp . The above property is obviously true for the 
present spin system, where E2 — i^min = but it should 
also remain valid for any realistic macroscopic system 
with non-degenerate ground state — consequence of the 
fact (jCl|) that N is exponentially larger than Ns- 

Indeed, let us consider as a more realistic example, 
the system of Ng spins 1/2 on a cubic lattice with ferro- 
magnetic nearest-neighbor exchange interaction and pe- 
riodic boundary conditions. The Hamiltonian is 7i = 
~JJ2nn ^™ ■ where Sm are vector spin operators, J 
is the (positive) exchange constant, and NN refers to the 
summation over the nearest neighbor {m, n}-pairs. In 
this case, the ground state has energy Ey^in = —jJNg. 
It is TVs-fold degenerate, because all states with total 
spin ^Ng have this energy. The excited states have 
a character of spin waves. Yet these spin waves are 
gapped by the finite linear size of the system {^^/N^), 

i.e. E2 - Smin ^ > As far as the ground 

state degeneracy is concerned, one can imagine, that any, 
even unrealistically small stray magnetic field would lift 
this degeneracy with the energy splitting that would still 



be exponentially larger than 



l-^minl 
N 



(It is also not diffi- 



cult to extend the result of this work to the theoretical 
case of degenerate ground state.) 



I begin with several general observations, which are 
applicable to both cases (i) and (ii). 

1) Let us denote as /i, the sum over the populations of 
those exceptional low-lying levels not describable by for- 
mulas (|18I20[) and exhibiting large average occupations 
(pk)- As mentioned in Section IIIDI there can exist only 
a relatively small number of these levels. Therefore, their 
energies should all be close to i?min- 

All other levels — in particular those surrounding the 
dominant peak in the density of states i^iE) — would 
have "regular" values of {pk) given by formula (^0)1 . 
Given the overall normalization (HI]), the total occupa- 
tion of regular levels, to be denoted as I2, is 



1-/1. 



(Dl) 



I can now use the energy condition and the fact 
that the occupations of the exceptional levels are peaked 
around -Bmim while the integrand -. — . fL^^^ — >■ is also 

peaked around Ec ~ [see Eq. ljClip and its derivation], 
to obtain 



{E„ 



Ea,v)Il — 



E, 



l-XE, 



-I2 = 0, 



(D2) 



which, in combination with (|Dip and the approximation 
((32| for A gives 



Em 



E„ 



E,: 



E^ 



(D3) 



(D4) 



APPENDIX D: CONDENSATION INTO THE 
LOWEST LEVEL OF MACROSCOPIC SYSTEM 

When condition ([28| is not fulfilled, volume Vk{pk) 
does not decay exponentially fast. The possible alter- 
natives are (i) that it changes slowly — decreases or in- 
creases, or (ii) it has a sharp Gaussian-like maximum 
at Pk = Pko somewhere between and 1, in which case 
(pk) ~ Pko- In the both cases, (pk) would be a significant 
fraction of 1. 

The first of the above cases might be realized in nu- 
merical studies, when N 1, but logA^ is not too large. 
In this case, several low-lying levels may exhibit large 
values of (pk)- 

I now show that in the macroscopic system, of the type 
described in Sectior HTEl it is the second of the above 
cases that is realized, and that the significant average 
occupation builds only for the lowest-lying level, while for 
all other levels, continue exhibiting occupations (pk) <C 1 
describable by formula ((20| . 

The subtlety of the present part is that it requires one 
to find a correction to a very accurate approximation 
(IM)) . In doing so, one needs to go back and forth be- 
tween the scale of jii'niinl and the scale of the level spacing 
between the lowest two levels. 



2) When the general formula ([31]) is used, one has to 
remember that manifold Mk exists and its volume Vk has 
non-zero value, only when < -Eav < -^max' where 

E'i„ and E'„^ are the minimum and the maximum value 
among all energies excluding Ek- This means that the 
argument of X[E] should also stay within the same limits, 
which, in turn imposes cutoff on Vk{pk) as a function of 
Pk, when the upper integration limit in Ea. (|3ip reaches 
E'^^^ or E[^.^. Beyond this cuttoff, Vk{pk) = 0. 

Here I am primarily concerned with the exceptional 
low- lying levels (and also exclude the case i?av < E2), 
which translates into constraint 



p _ {Ek - £^av)Pfc p 

1 ~Pk 

This results in the upper maximum value for pk '■ 



PkU 



E ~ E 

-Emax Eh 



(D5) 



(D6) 



One can establish a stronger effective cutoff for pk, 
beyond which Vk{pk) is guaranteed to exhibit sharp ex- 
ponential decay (for Ek < -Eav)- This cutoff corresponds 
to 



E,. 



{Ek — Eiiy)pk 

I- Pk 



0, 



(D7) 
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which translates into the cuttoff value 



I obtain after some manipulation: 



E.. 



PkC 



av 



(D8) 



The above cuttoff originates from the following argu- 
ment. When Vk{pk) given by ((3T|) decays slowly or in- 
creases at small p^, it happens, because the increasing 
value of integral in pip nearly compensates or outweighs 
the decreasing value of the preceeding logarithmic term. 
I note now, that X[E] becomes negative above E ~ 0. 
Therefore, once pk reaches the value given by Eq. (jD7|l . 
the integration extends into positive E and the integral 
starts decreasing. After that, it is certain that nothing 
any longer can slow the fast exponential decay — hence 
cutoff (IPSl). 



3) Cutoff (|D8p can now be used to show that as E'av 
decreases, the value of E\[Eiiy] cannot always stay below 
-Ernin- In principle, the estimate (jClip already indicates 
this, but it could have happened that it signified only Ex 
approaching -Bmin within 0{1/N) from below and never 
crossing it. 

If the latter possibility were to be realized, the deriva- 
tive of Vi{pi) would remain negative or zero for < 
Pi < 1 with rapid exponential drop above pic given by 
Eq. ljDSp . This would, in turn, imply that 



(Pi) < Ipic + Oil/N) 



1 E,. 



2 En 



Given inequality ([39]), (^2) < 1 for Ex < Emin- There- 
fore, the only "exceptional" level in this case would be 
the lowest one. That is, Ii = (pi), where {pi) is given 
by Eq.fUg)). which contradicts to Eq. (|D3ll . Thus the es- 
timate (IC11|) indicated correctly that Ex crossing E'min 
corresponds to E^^ — Ec ^ -Emin- 

4) When Ex > i?min, the volume Vi{pi) increases at 
small pi ; then it reaches a maximum at the value of pi to 
be denoted at pio, and, finally decays to zero at pi = 1. 

Below I locate piQ and show that it corresponds to a 
(5-function-like maximum: 



Vi{pi) exp 



(pi ~Pio)^ 



2cti 



where y^ctT ^ 1 , which implies 

{Pi) ~Pio- 



(DIO) 



(Dll) 



Requiring the derivative of the exponential power in 
Eq. ([3T|l with respect to pk to be equal to zero, I obtain 
for fc = 1: 



1 + A 



(Emin — E^^)p 



10 



1 -Pw 

Making in (jD12[) substitution 

m 



E ■ 



E,. 



1 - Pw 



E-Ex[Ey 



= 0. (D12) 



(D13) 



Ex 



E,: 



{Emin — Ei^^)p 



10 



1 



PlO 



E„ 



(DM) 



Therefore, according to the definition of Ec given by 
Eq.dSOD, 



E, 



{En 



Eav)pi 



As a result, 



Pw 



1 -Pio 



Ec — Eg, 
Ec - En, 



Ec- 



E, 



En 



(D15) 



(D16) 



The transition to the approximate value above follows 
from inequality 

The approximate value in (jPlSP is equal to the effec- 
tive cutoff value pic given by (|D8[) . where the derivative 
of Vk{pk) is supposed to be strongly negative. There 
is, however, no contradiction here. The accurate value of 
Pio given by the middle expression in Eq. (|D16[) is smaller 
than Plc^ and, as shown below, the width of the max- 
imum ^Joi is much smaller than the difference between 
Pic and piQ. 

The value of ui can now be obtained from 



0{l/N). (D9) 



1 



(P\ogVi{pi 



dpi 



(TV -3) dEx[E] 



dE 



E=Ec 

(D17) 

Condition — '^^J'i^^ ^ 1/A^ would then be sufS- 

E=Ec 

cient to justify cji ^ 1. 

By differentiating the accurate version of conditions 
(|D1|D2|1 and using formula ([31]) for Vi{pi) to calculate 
Ii = (pi), I was able to obtain 



dEx[E] 



dE 



2Ei 



Ns » 1, 



(D18) 



E=Ec 



which guarantees the validity of Eq. (jD9[) with pio ap- 
proximated by Eq. (|D16p — hence 



{Pi) 



E,, 



En 



(D19) 



Finally, I observe, that Eqs. ((D3)l and (ID19|) imply that 



{Pi) « h- 



(D20) 



The approximate values (|D3p and (|D19p for Ii and {pi) 
were obtained under the same assumption Ec ~ 0. In 
fact, Eq. (|D20[) also holds, when Ec is not neglected. 

Equation (jD20[) indicates that (pi) exhausts or almost 
exhausts the total occupation Ii of all exceptional levels, 
whose respective Hilbert space volumes depart signifi- 
cantly from formula 

One can further rule out the possibility "almost ex- 
hausts" . It would require that Ex [i?av] reaches at least 
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the second lowest energy level E2 ■ If this were to happen 
at ii^av — Ec2, then V2{p2) would aquire maximum at 



P20 



Ec2 — 



E, 



C2 



E„ 



(D21) 



for all Eav < Ec2- This maximum, sharp or not, would 
then lead to (pi) + (p2) becoming significantly greater 
than E^v/Eniin for most of the range E^v < £-02 in con- 
tradiction with Ea. (jD3|) . 

Therefore, E\ is always sufficently smaller than E2 to 
justify the use of formula ((20)) for all /c > 2. 



APPENDIX E: DERIVATION OF EQS. (|69||72|) 
FOR THE DENSITY MATRIX OF A 
SUBSYSTEM IN A MACROSCOPIC 
ENVIRONMENT 



It follows from inequality ([55)) that Ex[E] has the fol- 
lowing approximate property: 



and 



if Ex[E] < Enun, then E « Ep, (El) 



if S,nin <E<Ep, then Ex[E] « (E2) 



where Ep is the position of the sharp maximum of i^^E). 
It was also obtained in Appendix [Pl that 



Ex[E]<E2, yE<Ep. 



(E3) 



where E2 is the energy of the second lowest level. 

The key assumption of the present calculation is that 
the above conditions apply to both sets of energies {Ea} 
and {Ej,} with the appropriate insertion of subscripts a 
and b. The densities of states i^a{E) and i^biE) corre- 
sponding to these two sets are characterized by the re- 
spective peak positions: 



Epa = Esa, 
Esa 



^Pb 



Ni - 1 



(E4) 



(E5) 



Each energy set also has a respective minimum value 
and -Efjinin The outcome of the calculation now 
simply depends on whether -Eamin is smaller or larger 
than E^bmin- 
Case I: a > 2 
In this case: 



E„ 



E. 



Sa 



E 



Emin ; 



(E6) 



and 



Ebmin — Esniin + ^^Bmin = ^^min, (E7) 
which implies that -Eamin > ^-bmin- 



Given conditions ()E3[) and ([68)1 , 1 can further constrain 

E\t, < Eb2 < Eamin (E8) 

This condition applies to Ex^ [E] independently of its ar- 
gument. Therefore, 



Ek 



Ea: 



Ea 

Pa 



< Ean 



Thus, according to ()E1[) . 

Ea 



E,: 



Pa 



Epa — Esa- 



(E9) 



(ElO) 



One can now solve the pair of equations (|65IE10[) with 
respect to Ea and pa to obtain 



(Ell) 

jed by 
(E12) 



\ Eav — 


Ex^ 


Ea.v - 


I Ea 

Pa _ 


^Sa 


Ex^ 


Eav - 


I Ea 

Pa _ 


to find Ex^ 


Ea.v ' 


I Ea' 

Pa 


, one 



the inequality 



Eav < Epa,Epb, 



which is the consequence of the earlier assumption ()67p 
that |i?av| > Esa, Va. Equations (|E10|E12[) now require 
that Ea > 0. As a result. 



E, 



Ea 



< Eav < Epb, 



1 - Pq 

which, according to (|E2[) . implies that 

Ea 



Ex 



Ea 



1 - Pa 

Therefore, according (|66)l and ()E7I 



= Ebn 



Ex. 



-C/av H 

Pa . 



(E13) 



(E14) 



(E15) 



The substitution of Eq. (|E15P into Eq. ijEllj) yields 

1 Eav ~ £"111111 l-c-ya\ 

Po^ = JT-f^ 7? ■ (E16) 

^Vl rjSa — Aniii 

Given inequality ()67p . the leading order approximation 
of (|El6)) is 



1 A Ea 



(E17) 



Case II: a = 1. 
Now 



-^amiii — -^Siiiiii -£i?min — E-^ 



(E18) 
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and 



Ebmin — E, 



S2 



Emin i 



(E19) 



Given the condition |i?smin| ^ |-E'av|, l-^^Bminl, the lead- 
ing order approximation of (|E25p is 



implying Samin < -Bb,ni„. 

One can obtain using condition (|68p . that -Bftmin — 

Eaiaiu = Es2—Esmiu > EE2 — EEmin — Ea2~Eamm- The 

resulting inequality together with condition (|E3p implies 
that 



E, 



Pi 



E 



Smin 



1 - 



E, 



E 



En\\i\ 



(E26) 



E\^ < Ea2 < -Efcmin- (E20) 

Inequality (|E20p applies to any argument of E\ [E] in- 

Ea 
Pi 



eluding i?av + Using this fact together with Eg. ([66|) . 



I obtain 



Ex, 



E,. 



Ea 



I -Pi 



<E, 



6min • 



Thus, according to (|Eip . 

Ea 



E,: 



I -Pi 



Epb — — 



Esn 



Ni - 1 



(E21) 



(E22) 



Given inequality (|E12p . Equation (IE22P implies that 
Ea < 0. As a result, 



-Eav H — < Eav < Epa, 

Pi 



and thus, according to (|E2 



E^ 



E, 



Ea 
Pi 



E„ 



(E23) 



(E24) 



After Eq. (|E24p is substituted into ([65]) and the result- 
ing equation is solved jointly with (|E22p with respect to 
pi and Ea , I obtain 



Pi 



-Eon 



7Vi-l 



Afi-1 



Err 



(E25) 



It is clear from the derivation that formulas (jE16|E25P 
describe a broader range of mostly abstract possibilities 
beyond the physical assumptions used so far. 

One can, in particular, consider relaxing the condition 
\Esa\ ^ l-E-avl, l^'Bminl- In this casc, however, on needs 
to be concerned with the validity of conditions (|E11 IE2p 
for the energy set {Ei,}. For example, if the spectrum of 
{Esa} consists of only a few far separated levels, then 
Vb{E) loses the single peak structure and hence condi- 
tions (|E11 IE2p . However, when iVi ^ 1, and the spec- 
trum of energies {Esa} has a sharply-peaked Gaussian- 
like density of states, then conditions (|Eli IE2p are recov- 
ered. 

In the latter case, one also needs to examine what hap- 
pens, when i?av falls between Epa and i?pb. In this case, 
the system of equations (I65[ I66p with approximations 
(IE11IE2P generates two more solutions, in addition to the 
one already found. These two solutions correspond to the 
possibility that both E^^ + ^ and i?av — become 
approximately equal to Epa and Ep^ respectively. One 
solution would then correspond to positive values and 
the other one to the negative values of Aa, Af,. They are 
presumably a maximum and a minimum located close 
to each other in the plane of variables {EA,Pa)- The 
additional maximum should, presumably give a smaller 
Hilbert space maximum, in comparison to the one found 
earlier, because it has a larger value of pa, and there- 
fore stronger suppressed by term containing log(l — Pa) 
in Eq.([5D|. 
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